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Abstract 

The notion of stability can be generalised to point processes by defining the 
scaling operation in a randomised way: scaling a configuration by t corresponds to 
letting such a configuration evolve according to a Markov branching particle system 
for -logf time. We prove that these are the only stochastic operations satisfying 
basic associativity and distributivity properties and we thus introduce the notion 
of branching-stable point processes. We characterise stable distributions with re¬ 
spect to local branching as thinning-stable point processes with multiplicities given 
by the quasi-stationary (or Yaglorn) distribution of the branching process under 
consideration. Finally we extend branching-stability to random variables with the 
help of continuous branching (CB) processes, and we show that, at least in some 
frameworks, ^-stable integer random variables are exactly Cox (doubly stochastic 
Poisson) random variables driven by corresponding CB-stable continuous random 
variables. 
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1 Introduction 

The concept of stability is central in Probability theory: it inevitably arises in various 
limit theorems involving scaled sums of random elements. Recall that a random vector 
£ (more generally, a random element in a Banach space) is called strictly a-stable or 
StaS, if 

t 1/a e + (1 - t) 1/a £" = £ for all t G [0,1], (1) 
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where £' and £" are independent copies of £ and = denotes equality in distribution. 
When a limiting distribution for the sum of n independent vectors scaled by n l ^ a exist, 
it must be StaS, since one can divide the sum into first tn and the last (1 — t)n terms 
which, in turn, also converge to the same law. This simple observation gives rise to 
the defining identity (jTj) . It is remarkable that the same argument applies to random 
elements in much more general spaces where two abstract operations are defined: a sum 
and a scaling by positive numbers which should satisfy mild associativity, distributivity 
and continuity conditions, i.e. in a cone, see [6]. For instance, a random measure £ on 
a general complete separable metric space is called strictly a-stable if identity (JTJ) is 
satisfied, where the summation of measures and their multiplication by a number are 
understood as the corresponding arithmetic operations on the values of these measures 
on every measurable set. Stable measures are the only class of measures which arise as 
a weak limit of scaled sums of random measures. 

Since the notion of stability relies on multiplication of a random element by a number 
between 0 and 1, integer valued random variables cannot be StaS. Therefore Steutel 
and van Harn in their pioneering work [19] defined a stochastic operation of discrete 
multiplication on positive integer random variables and characterised the corresponding 
discrete a-stable random variables. In a more general context, the discrete multiplication 
corresponds to the thinning operation on point processes when a positive integer random 
variable is regarded as a trivial point process on a phase space consisting of one point (so 
it is just the multiplicity of this point). This observation leads to the notion of thinning 
stable or discrete a-stablc point processes (notation: DaS) as the processes <f> which 
satisfy 

+ (1 -t) 1/a o&' = $ for all! e [0,1], (2) 

when multiplication by a f G [0,1] is replaced by the operation to of independent thinning 
of their points with the retention probability t. The DaS point processes are exactly the 
processes appearing as a limit in the superposition-thinning schemes (see [H Ch. 8.3]) 
and their full characterisation was given in |T]. 

In its turn, a thinning could be thought of as a particular case of a branching op¬ 
eration where a point either survives with probability t or is removed with the com¬ 
plimentary probability. This observation leads to a new notion of discrete stability for 
point processes by considering a more general branching operation based on a subcriti- 
cal Markov branching process (Y t ) t> 0 with generator semigroup T = (F t ) t > 0 , satisfying 
Y 0 — 1. Following Steutel and Van Harn [21j who considered the case of integer-valued 
random variables, we denote this operation by o T . In this setting when a point process is 
“multiplied” by a real number t G (0,1], every point is replaced by a collection of points 
located in the same position of their progenitor. The number of points in the collection 
is a random variable distributed as Y_i ogf . This operation preserves distributivity and 
associativity with respect to superposition and generalises the thinning operation. In 
Section [3] we study stable point processes with respect to this branching operation o T 
calling them ^-stable point processes. We show that ^-stable point processes are es¬ 
sentially DaS processes with multiplicities which follow the limit distribution Yqo of the 
branching process Y t conditional on its survival (Yaglom distribution) and we deduce 
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their further properties. 

In a broader context, given an abstract associative and distributive stochastic oper¬ 
ation • on point processes, a process $ is stable with respect to • if and only if 

Vn G N 3c n G [0,1] : <L = c n • ($ (1) + ... + <L (n) ), 

where ; $( n ) are independent copies of <h. In such a context stable point processes 

arise inevitably in various limiting schemes similar to the central limit theorem involving 
superposition of point processes. In Section 0] we study and characterise this class of 
stochastic operations. We prove that a stochastic operation on point processes satisfies 
associativity and distributivity if and only if it presents a branching structure: “multiply¬ 
ing” a point process by t is equivalent to let the process evolve for time — log t according 
to some general Markov branching process which may include diffusion or jumping of the 
points. We characterise branching-stable (i.e. stable with respect to •) point processes 
for some specific choices of •, pointing out possible ways to obtain a characterisation 
for general branching operations. In order to do so we introduce a stochastic operation 
for continuous frameworks based on continuous-state branching Markov processes and 
conjecture that branching stability of point processes and continuous-branching stability 
of random measures should be related in general: branching-stable point processes are 
Cox processes driven by branching-stable random measures. 

2 Preliminaries 

In this section we fix the notation and provide the necessary facts about branching 
processes and point processes that we will use. We then address the notion of dis¬ 
crete stability for random variables and point processes that we generalise in subsequent 
sections. 

2.1 Branching processes refresher 

Here we present some results from 0 ch.m], mi ch.v], m and eh about continuous 
branching processes that we will need. Let (W) s > 0 be a Z + -valued continuous-time 
Markov branching process with Y 0 — 1 almost surely, where Z + denotes the set of non¬ 
negative integers. Markov branching processes describe the evolution of the total size of 
a collection of particles undergoing the following dynamic: each particle, independently 
of the others, lives for an exponential time (with fixed parameter) and then it branches , 
meaning that it is replaced by a random number of offspring particles (according to 
a fixed probability distribution), which then start to evolve independently. Such a 
branching process is governed by a family of probability generating functions (p.g.f.’s) 
T = (F s ) s >o, where F s is the p.g.f. of the integer-valued random variable Y s for every 
s > 0. It is sufficient for us here to consider the domain of F s to be [0,1]. It is well 
known that the family F is a composition semigroup: 

F s+t (-) = F s (F t (-)) Ws,t> 0. (Cl) 
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Conversely, by writing relation (jC 1 j) explicitly in terms of a power series, it is straight¬ 
forward to see that the system of p.g.f.’s corresponding to an non-negative integer valued 
random variables describes a system of particles branching independently with the same 
offspring distribution, so that Y s is the number of particles at time s, see also HD 
Ch.V.5], 

We require that the branching process is subcritical , i.e. E[Y S ] = F s '(l) < 1 for s > 0. 
Rescaling, if necessary, the time by a constant factor, we may assume that 


E|yj = y<i) = e-. 

Finally we require the following two regularity conditions to hold: 

0 < 2 < 1 , 


limF s (^) = F 0 (z) = z 
lim FJz) = 1 


0 < z < 1. 


(C2) 

(C3) 

(C4) 


(1C-ill implies that the process starts with a single particle Y 0 = 1 and (1C4|) is a conse¬ 
quence of the subcriticality meaning that eventually Y s = 0. 

A rationale behind requiring (1C2|) . (IC3[) and (1C4D will become clear later, see Re¬ 
mark [H Identities (1CID and (IC3D imply the continuity and differentiability of F s (z) 
with respect to s, see, e.g., [21 Sec.III.3], and thus one can define the generator of the 
semigroup T 


U ^ = T S 


FJz) 


0< z<l. 


s =0 


The function U(-) is continuous and it can be used to dehne the A-function relative to 
the branching process 


A(z): — exp 


dx 

Ufx) J 


0 < z < 1, 


(3) 


which is a continuous strictly decreasing function with A(0) = 1 and A(l) = 0, see, 
e.g., [21 Sec.III. 8]. From (jCljl it follows that U(Fjz)) = U(z)Fjz) and therefore 


A(F s (z)) = e s A(z) s > 0, 0 < z < 1. 


(4) 


Definition 1. Let (W) s > 0 and T = (R s )s> o be as above. The limiting conditional 
distribution (or Yaglom distribution) of Y s is the weak limit of the distributions of 
(W|Y S > 0) when s — > +cx). We denote by Y 0 0 the corresponding random variable and 
by B(-) its p.g.f., called the B-function of Y s . 

The B-function of Y s is given by 

B{z): = l-A(z)= lim 0 < z < 1. (5) 

From © and ([5]) it follows that 

B(Fjz)) = 1 — e~ s + e~ s B(z), s > 0, 0 < z < 1. (6) 
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Both A and B are continuous, strictly monotone, and surjective functions from [0,1] 
to [0,1], thus the inverse functions A _1 and B~ 1 exist and have the same properties. 
Moreover, using fl5]) we obtain 

e~ s = 1. 

s=0 

At the same time 


a 


-rA(FM) 


.9=n 


a 

ds 


d_ 

ds 


A(F,( 0))=V(0)h 
■0 ds 


F,( 0) = A'( 0)h 

o ds 


P{Y(s) = 0}. 

s=0 


Since Y s is a continuous Markov branching process, every particle branches after expo¬ 
nentially distributed time with a non-null probability to die out and it follows that 


d_ 

ds 



0 } > 0 


implying also that 


A'(0) = 


d 


.ds 


s =o 


-l 


P— 0} G (0, +oo). 


(7) 


The simplest, but important for the sequel example is provided by a pure-death 
process. 

Example 1. Let (Y s ) s>0 be a continuous-time pure-death process starting with one indi¬ 
vidual 

Y = 1 1 if s < r ’ 

| 0 if s > t, 


where r is an exponential random variable with parameter 1. The composition semi¬ 
group T = (F s ) s>0 driving such a process is 


F s {z) = 1 - e~ s + e~ s z 0 < z < 1. (8) 

Clearly T = (Ps) s>0 satisfies (1C II) iK'lli . The generator U(z) and the A and 5-functions 
defined above are 


U(z) = A(z) = 1 — z, B(z) — z, 0 < z < 1. (9) 

Another example is the birth and death process. 

Example 2. Given two positive parameters A and //, assume that each particle disappears 
from the system at rate // or it is replaced with two particles at rate A independently of 
the others. The total number of particles at each time can either grow or diminish by 
one, hence this process is also called the linear branching. Its generator is given by 

U(z) — n — (A + n)z + A z 2 . (10) 
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The process is subcritical whenever p > A and in order to satisfy (1C2|i one needs to scale 
the time so that /i = A + 1. This defines a one-parametric family of semigroups 


F,{z) 


e S (1 — z) 

1 + A(1 — e _s )(l — z) ’ 


see pj p. 109]. 
given by 


Conditions (1C1I) , (1C3I) and (1C4[) also hold and the functions A and B are 


A{z) 


(A + 1)(1 - z) 

1 + A(1 — z) ’ 


B{z) 


z 

1 + A(1 - z) ’ 


for z in [0,1]. Thus B describes the p.g.f. of a (shifted) Geometric distribution with 
parameter (1 + A) -1 . 


2.2 Point processes refresher 

We now pass to the necessary definitions related to point processes. The details can be 
found, for instance, in a. a and m- A random measure on a phase space X which 
we assume to be a locally compact second countable ffausdorff space, is a measurable 
mapping £ from some probability space (P, T. P) into the measurable space (At, B( At)), 
where At denote the set of all Radon measures on the Borel a-algebra B(X) of subsets 
of X and B( At) is the minimal cr-algebra that makes the mappings p K > p(B), p £ At 
measurable for all B £ B(X). 

The distribution of a random measure is characterised by the Laplace functional 
L^[u] which is defined for the class BM + (A) of non-negative bounded functions u with 
bounded support by means of 

L^[u\ = Eexpj — (w,£)}, m £ BM + (T). (11) 

Here and below (u,p) stands for the integral f u(x ) p(dx ) over the the whole X unless 
specified otherwise. 

A point process (p.p.) $ is a random counting measure, i.e. a random measure that 
with probability one takes values in the set J\f of all boundedly finite counting measures 
on B(X). The corresponding a-algebra B(Af) is the restriction of B(M) onto J\f. The 
support of p £ J\f is the set supp(^:={x £ X : <^({a:}) > 0}. A point process is 

called simple (or without multiple points) if P<|) _1 {(^ £ J\f : (^({a:}) < 1 \/x £ X} — 1. 
The distribution of a point process $ can be characterised by the probability generating 
functional (p.g.fl.) G$[h] defined for functions h such that 0 < h(x) < 1 for all x E X 
and such that the set {x £ X : h(x) ^ 1} i s compact. We denote the class of such 
functions by V(X). Then 

G$[h] = L$[— log h] = E exp {(log h, <h) j, h £ V(X). 

For a simple p.p. $, this expression simplifies to 

G*[h]= E Yl K x i)i heV(X). 

iCiGsupp $ 
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A Poisson point process with intensity measure A is the p.p. II having the p.g.fl. 


Gu[h] = exp{ — (1 — h, A)}, hEV(X). 


It is characterised by the following property: given a family of disjoint sets B j G B(X), 
i = 1,..., n, the counts II(5i),..., II( ) are mutually independent Poisson Po(A(5j)) 
distributed random variables for i — 1 , ,n. 

Given a random measure £, a Cox process with parameter measure f is the point 
process with the p.g.fl. 


G*[/i] = Eexp{-<l-./i,f)}, heV{X). 


( 12 ) 


It is called doubly-stochastic, since it can be constructed by first taking a realisation £(cu) 
of the parameter measure and then taking a realisation of a Poisson p.p. with intensity 
measure £(cu). 

Consider a family of point processes ('Pj / ) yg y on X indexed by the elements of a locally 
compact and second countable Hausdorff space y which may or may not be X itself. 
Such a family is called a measurable family if P y (A) : = P(^ y G A) is a £>(jb)- measurable 
function of y for all A G B{M). 

Given a point process S on y and a measurable family of point processes {^y)yey 
on X, the cluster process is the following random measure: 



(13) 


The p.p. S is then called the center process and y G y are called the component 
processes or clusters. The commonest model is when the clusters in (I13p are independent 
for different y\ G suppS given a realisation of S. In this case, if G~\h] is the p.g.fl. of 
the center process and Gq,[h\y\ are the p.g.fl.’s of y G y, then the p.g.fl. of the 
corresponding cluster process (Q~3j) is given by the composition 


G^[h] = G s [G*[/i|-]]- 


(14) 


2.3 Stability for discrete random variables 

Let A" be a Z + -valued random variable. As in [T2j, we define an operation of discrete 
multiplication o by a number t G [0,1] 



i =1 


where {Zb)} igN are independent and identically distributed (i.i.d.) random variables 
with Bernoulli distribution Bin(l,t). A random variable X (or its distribution) is then 
called discrete a-stable (notation: DckS) if 


t 1/a o X' + (1 - t) 1/a o X" = X for all t G [0,1], 


(16) 
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where X', X" are independent distributional copies of X. 

Letting each point i evolve as a pure-death process independently of the others 
(see Example [Tj) , after time — logt, t G (0,1], the number of surviving points will be 
distributed as (IEjI) . So alternatively, 


1=1 

Replacing here the pure-death with a general branching process allowed the authors 
of [21] to define a more general branching operation and the corresponding F-stable 
non-negative integer random variables as follows. Let be a sequence of i.i.d. 

continuous-time Markov branching processes driven by a semigroup T = (F s ) s >o satis¬ 
fying the conditions (C1)-(C4) in the previous section. Given t G (0,1] and a Z + -valued 
random variable A" (independent of {E^}ieN) define 

torX: = f2 Y - V ( ir > 

2=1 

The notion of ^-'-stability is then defined in an analogous way to discrete stability: 

Definition 2. A Z + -valued random variable X (or its distribution) is called T-stable 
with exponent a if 


t 1/a ojr X' + (1 - t) l/a o T X" = X Vt G [0,1], (18) 

where X' and X" are independent copies of X. 

In terms of the p.g.f. Gx(z) of X, (fT8|i is equivalent to 

Gx(z) = Gx(F_ a -i\ 0 g t (z)) ■ G j jf(E_ Q ,-i l0 g(i_t)(z)) 0 < z < 1. 


Let G tOT x{z) denote the p.g.f. of t o r X. By independence of {hW(-)} ieN and X, 
(HID is equivalent to 


Gtojrx(z) — Gx(F^i ogt (z)) 0<z<l. (19) 

ft is easy to verify that (1CID and dT9]l make the branching operation o T associative, 
commutative and distributive with respect to the sum of random variables, i.e. for all 
t,ti, t 2 G [0,1] and X independent of X' 

tl 0 F (h °T A") = (tlt 2 ) O jr X = t 2 Ojr (tl Ojr A"), 

t o T (A + X') =tOfX + tOf X'. 
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( 20 ) 

( 21 ) 



Remark 1. As shown in [2(1, Section V.8, equations (8.6)-(8.8)], conditions (1C2|) . (1C3|) 
and (IC4f) guarantee that o T has some “multiplication-like” properties. In particular 
(iC3l) and (1C4[) imply respectively that lim^] t.o T X = lo T X = X and lim,y 0 1ojr A = 0. 
Furthermore, (1C2[) implies that, in case the expectation of A" is finite, EftojpA"] — t EA. 

The following theorem gives a characterisation of .F-stable distributions on Z + , see 
[21, Theorem 7.1] and [2D;, Theorem V.8.6]: 

Theorem 1. Let X be a r L + -valued random variable and Gx{z) its p.g.f., then X is 
T-stable with exponent a if and only if 0 < a < 1 and 

G x (z ) = exp { - cA(z) a } 0 < z < 1, 

where A is the A-function (J3[) associated to the branching process driven by the semigroup 
T and c > 0. In particular, X is DoS if and only if 

G x (z ) = exp { - c(l - z) a } 0<z<l 

for some 0 < a < 1 and c > 0, see III). Theorem 3.2]. 

2.4 Thinning stable point processes 

It was noted in Section 12.31 that the operation of discrete multiplication of an integer 
random variable to with t 6 [0,1] may be thought of as an independent thinning when 
the random variable is represented as a collection of points and each point is retained 
with probability t and removed with the complementary probability. Thus the thinning 
operation generalises the discrete multiplication to general point processes. The corre¬ 
sponding thinning-stable or discrete a-stable point processes (notation: DaS) satisfy ([2]) 
and are exactly the ones which appear as the limit in thinning-superposition schemes, 
see HU Ch.8.3]. The full characterisation of these processes is given in [7], Thinning 
stable processes exist only for a G (0, 1], and the case a = 1 corresponds to the Poisson 
processes. 

To be more specific, we need some further definitions. First we need a way to 
consistently normalize both finite and infinite measures. Let be a fixed 

countable base of the topology on X that consists of relatively compact sets. Append 
B 0 = X to this base. For each non-null // G A4 consider the sequence of its values 
(/i(-B 0 ), p-(Bi), /a(B 2 ),...) possibly starting with infinity, but otherwise finite. Let i(p) 
be the smallest non-negative integer i for which 0 < p(Bi) < oo, in particular, i(p) = 0 
if p is a finite measure. Define 

§ = {/j, G M : ia{B i{fl) ) = 1}. 

It can be shown (see 0) that S is B(M )-measurable and that Sfl{/r : p(X) < oo} = Mi 
is the family of all probability measures on X. Furthermore, every /i G Xi \ {0} can be 
uniquely associated with the pair (/t, p(Bi^)) G § x M + , where ft is defined as p > 
and n = /i(5j( M ))/i is the polar representation of (i. 
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A locally finite random measure £ is called strictly stable with exponent a or StaS if 
it satisfies identity ([Tj) . It is deterministic in the case a = 1 and in the case a e (0,1) 
its Laplace functional is given by 

L^[h]=ex pi- f (1 - e~ M )A(dp)}, he BM+(A), (22) 

^ Jm\{ 0} J 

where A is a Levy measure, i.e. a Radon measure on Ad \ {0} such that 

f (i - e -(M)A(d//) < oo (23) 

Jm\{ 0} 

for all h e BM + (A). Such A is homogeneous of order —a, i.e. A (tA) = t~ a A(A) for all 
measurable A C M \ {0} and t > 0, see ]7J Th. 2], 

Introduce a spectral measure a supported by § by setting 

a (A) = T(1 — a) A({tp : p e A, t > 1}) 

for all measurable A C S, where T is the Euler’s Gamma-function. Integrating out the 
radial component in (1221) leads to the following alternative representation [7, Th. 3]: 


L^[u] = expi — / (u, p) a a(dp) \ , u e BM + (A) 


(24) 


for some spectral measure a supported by S which satisfies 


p(B) a a(dp) < oo 


(25) 


for all relatively compact subsets B of X. The latter is a consequence of (1231) and 
representation (1241) is unique. 

The importance of StaS random measures is explained by the fact that any DckS 
point process <f> is exactly a Cox processes driven by a StaS parameter measure £: its 
p.g.fl. has the form 

G$[h\ = L^[l — h] = expj— J(1 — h, /i)"cr(dp)|, h e V(A) (26) 

for some locally finite spectral measure a on S that satisfies (l25]h see [7) Th. 15 and 
Cor. 16]. 

In the case when a charges only probability measures M l5 the corresponding DcxS 
p.p.’s are cluster processes. Recall that a positive integer random variable 7] has Sibuya 
Sib(a() distribution with parameter a, if its p.g.f. is given by 

E ^ = 1 — (1 — z) a , 2 6(0,1], 

It corresponds to the number of trials to get the first success in a series of Bernoulli 
trials with probability of success in the A;th trial being a/k. 
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Definition 3. (See [7, Def.23]) Let p be a probability measure on X. A point process 
T on X defined by the p.g.fl. 

G r [h] = GrU h ] = 1 - (1 - h,p) a , h G V(X), (27) 

is called a Sibuya point process with exponent a and parameter measure p. Its distri¬ 
bution is denoted by Sib (a,p). 

A Sibuya process T ~ Sib(o,/r) is a.s. finite, the total number of its points T(X) 
follows Sib(a) distribution and, given the total number of points, these points are inde¬ 
pendently identically distributed in X according to p. 

Theorem 2 (Th. 24 0)- A DoS point process with a spectral measure a supported by 
Mi can be represented as a cluster process with Poisson centre process on Mx driven by 
intensity measure a and component processes being Sibuya processes Sib (a,p), p G 
Its p.g.fl. is given by 

G<s>[h] = exp| [ (G T ^[h} - 1) a(dp)\, h G V(X), 

VM, J 

with [h] as in (1271) . 

3 ^-stability for point processes 

We have seen in the previous section that the discrete multiplication operation on integer 
random variables generalises to the thinning operation on points processes. In a similar 
fashion, we can extend the branching operation o T to point processes too. 

Let {y s }. s >o be a continuous-time Markov branching process driven by a semigroup 
T = (F s ) s > o satisfying conditions (ICl|) - (|C4j) . Intuitively, given a point process <3> and 
f G (0,1], t ojr $ is the cluster point process obtained from $ by replacing every point 
with Y_ i ogt points located in the same position (using an independent copy of Y_ j ogt for 
each point including the ones in the same position). In this sense, the resulting process 
is a cluster process. To proceed formally, we first define its component processes. 

Definition 4. Given a Z + -valued random variable Z and x G X. we denote by Z x the 
point process having Z points in x and no points in A\{a;}. Equivalently Z x is the point 
process with p.g.fl. Gz x [h] = F(h(x)) for each h G V(X), where F(z) is the p.g.f. of Z. 

We can now define the operation o T for point processes. 

Definition 5. Let $ be a point process and {W}s>o be a continuous-time Markov 
branching process driven by a semigroup T = ( F s ) s > 0 satisfying conditions (1C ID - (IC4f) . 
For each t G (0,1], t ojr is the (independent) cluster point process with center process 
$ and clusters {(bAiog^, x G T}. 
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Equivalently t o j <f> can be defined as the point process with p.g.fl. given by 

G t o T <&[h} = G r $[F_i ogt (/i)], 

where is the p.g.fl. of <h. Note that does not need to be simple (i.e. it can have 

multiple points), even if <f> is. We define the ^-stability for point processes as follows. 

Definition 6. A p.p. $ is Testable with exponent a (a-stable with respect to oj P ) if 

fl/ a & + (1 - t) 1/a o T $ VtG (0,1], (28) 

where <f>' and <f>" are independent copies of <f>. 

Equivalently, d28|) can be rewritten in terms of p.g.fl.’s as follows: 

G*[h] = G*[A_ logt/a (h)] Gz[F_ log{1 _ t)/a (h)\ Vt G (0,1], Vh G V(X). 


Remark 2. The branching operation o j induced by the pure-death process of Example [[] 
corresponds to the thinning operation. Therefore DoS point processes can be seen as a 
special case of .A-stable point processes. 

An .A-stable point process <h is necessarily infinitely divisible. Indeed, iterating (128|) 
m — 1 times we obtain 

m~ 1/a o T + ... + m~ l/a op = $, (29) 

where $d),... ; $( m ) are independent copies of d>. 

A characterisation of .A-stable point processes is given in the following theorem which 
generalises [2, Tli.15] and which proof we largely follow here. 

Theorem 3. A functional (?$[•] is the p.g.fl. of an T-stable point process $ with expo¬ 
nent of stability a if and only if 0 < a < 1 and there exists a StaS random measure f 
such that 

G*[h] = [ A(h )] = Lf: [1 - B(h )] Vh G V, (30) 

where A(z) and B(z) are the A-function and B-function of the branching process driven 

by T. 

Proof. Sufficiency: Suppose (I5UD holds. As it was shown in [7], StoS random measures 
exist only for 0 < a < 1, a = 1 corresponding to non-random measures. Next, by CO 
and (fT2j) . L g [l — h\ as a functional of h is the p.g.fl. of a Cox point process with intensity 
£ and B(z) is the p.g.f. of the limiting conditional distribution of the branching process 
driven by T. Therefore, by (JHJ), the functional G$[h] — L^[l — B(h)\ is the p.g.fl. of 
a cluster process, say <h. We need to prove that <f> is .A-stable with exponent a. Given 
t £ (0,1] and h G V(X) it holds that 


log t/a (h)] G§ [F_ log(l -t)/a (A)] = 


= L, 


A(A_ \ og t/a(h)) L(: A(A_i og (i_q/ Q (h)) 


* L^Aa)} Ltia-t^Aw]. 
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Since £ is StaS, it satisfies (H]) and thus 

Ld tU ° A W\ L^l-t^AW] = L t {A(h)}. 

Therefore 

G,[F. 

log t/a. (h)] G$ [F_ log(l -t)/a (h)} = G*[h\ 
for any h in V(X), meaning that $ is J-'-stable with exponent a. 

Necessity: Suppose that $ is Testable with exponent a. Writing (129j) for the values of the 
measures on a particular compact set B e B(X), we see that $(-8) is an J-’-stable random 
variable with exponent a. Thus by Theorem |Tj we have 0 < a < 1. Now we are going to 
prove that G$[A _1 (m)], as a functional of u, is the Laplace functional of a StaS random 
measure. While a Laplace functional should be defined on all (bounded) functions with 
compact support, the expression G<$,[A~ l (u)\ is well defined just for functions with values 
on [0,1] because A _1 : [0,1] —> [0,1]. To overcome this difficulty we employ (1291) which 
can be written as 


G$[h] = ( G4F a - Hogm (h)}) m Vh e V(X), 

and define 

L[u\ = (G^[F a - Hogm (A-\u))]y i i (G*[A-\m- l / a u)]y i u e BM + (fL), (31) 

for any m > 1 such that m~ 1 ^ a u < 1 . Note that the right-hand side of (19T1) does not 
depend on m. Moreover, given < 1, the function A^ l {m^ l / a u) does take values 

in [0,1] and equals 1 outside of a compact set, implying that A~ l {m^ l ^ a u) G V(X). 
Therefore L[u] in (1911) is well-defined. Since (19Tj) holds for all m, it is possible to pass 
to the limit as m —> oo to see that 

L[u\ = exp { — lim m(l — G$[A _ 1 (m _ 1 /“u)])) (32) 

l m—>oo ) 

We now need the following fact: 

lim m(l - G^[A-\m- 1 /a u)}) = lim m(l - G*[e^' 1 > (0 > m-1/att ]). (33) 

ra—>• oo m—>• oo 

Since A _1 is continuous, strictly decreasing and differentiable in 0 with t4 _1 (0) = 1 
and (t4 _1 ) / (0) < 0 (see Section 12TTT) . it follows that for any constant e > 0 there exists 
M(e, u) > 0 such that 

A-\ m - 1/a u(l + e)) < e {A- l Y(o)m-V"u < _ e )) V m > M(e, u ). ( 34 ) 

From (1921) . QMD and the monotonicity of we can deduce 

L[( 1 + e)u] < exp {— lim m(l — b'(o)m 1/Qtl ]) j < £[(i _ e ) M ] . ( 35 ) 

l m—>00 ) 
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Note that L is continuous because of (l3Tf and the continuity of G®. Therefore taking 
the limit for £ going to 0 in (1351) and using (1321) we obtain (1331) . 

From (1371) and Schoenberg theorem [3] Theorem 3.2.2] it follows that L[u] is positive 
definite if lim rn _ ) . 00 m(l — G$[l — is negative definite, i.e. by (1331) if 

n 

V dCj lim mil - < 0, (36) 

' m—>oo 

i,j =1 

for all n > 2, for any Uq,... ,u n E BM + (5) and for any ci,..., c n with Yh c i — 0. If we set 
Vi = e^ A b(o)m 1/a ui, (-} ien (| 3 ^|) j s equivalent to Y^7j=i c * c i lim m ^oo G$\viVj] > 0 , which 
follows from the positive definiteness of Thus, by the Bochner theorem [3] Theorem 
4.2.9], the function L[J2i =1 tihi] of t\,... ,tk > 0 is the Laplace transform of a random 
vector. Moreover L[0] = 1, where 0 is the null function on X. Finally from (1311) and the 
continuity of the p.g.fl. Gq, it follows that given {/ n }neN C BM + (5), /„ | / 6 BM + (5) 
we have L[f n \ —> L[f] as n —> oo. Therefore we can use Theorem 9.4.II in [5] to obtain 
that L is the Laplace functional of a random measure £. 

In order to prove that £ is StaS, let u £ BM + (5) and take an integer m > (sup u) a 
and denote by u = m~ x ' a u < 1. By CD, for any given t £ ( 0 , 1 ] we have 

h\u\ = G™{A-\u)} ^ G™[F„ logt/a (A-\u))] G™[F_ log{1 _ t)/a (A-\u))] i 

G^[A-\t l > a u)] G%[A-\(1 - t)V a u)] = L^[t 1/a u\ Ls[( 1 - t)V a u], 

which implies that £ is StaS. □ 

Corollary 1. A p.p. $ on X is F-stable with exponent a if and only if it is a cluster 
process with DaS centre process T on X and component processes {(Too)®, x £ Xj (see 
Definitions [7] and ]). 

Proof. From Theorem [3] and ([5]) it follows that $ is Testable if and only if its p.g.fl. 
satisfies G$[h] = L^[ 1 — 5(h)], where 5(-) is the p.g.f. of Y 0Q , and £ is a StaS random 
measure. By (ITilj) there is a DaS point process T with Gq,[h] — L^[l — h]. We obtain 
that Gq>[h] = Gq, [5(h)]. The result follows from the form (ITT]) of the p.g.fl. of a cluster 
process. □ 

Corollary Q] clarifies the relationship between 5-stable and DaS point processes: 5- 
stable p.p.’s are an extension of DaS p.p.’s, where every point is given an additional 
multiplicity according to independent copies of (the latter is fixed by 5). Note that 
when the branching operation is thinning, the random variable Y is identically 1 (that 
stems from ([9])) and the 5-stable p.p. is the DaS centre process itself. 

Corollary 2. A p.p. <F is 5 -stable with exponent 0 < a < 1 if and only if its p.g.fl. 
can be written as 

G$[tt] = exp{ - [(l - B(u),n) a cr(dfj.)}, (37) 

Js 

where a is a locally finite spectral measure on § satisfying ([25]) . 


14 


Proof. If $ is an Testable point process with stability exponent a, then by Theorem [3] 
there exist a StaS random measure £ such that 

G*[h] = Lt[A{h)] heV{X). 

Thus ([37]) follows from spectral representation (T24l) . Conversely, if we have a locally finite 
spectral measure a on § satisfying (123]) and a £ ( 0 , 1 ], then er is the spectral measure 
of a StaS random measure £, whose Laplace functional is given by (122]) . Therefore (1371) 
can be written as 

G*[h] = Lt[l~B(h)], 

which, by Theorem [3] implies the ^-'-stability of <f>. □ 

We also get the following generalisation of Theorem [2J 

Theorem 4. An T-stable point process with a spectral measure a supported only by 
the set M x of probability measures can be represented as a cluster process with centre 
process being a Poisson process on Mi driven by the spectral measure a and daughter 
processes having p.g.fl. Gr(n)[B(h)}, where T(/i) areSib(a,/i) distributed point processes 
and B(-) is the B-function of the branching process driven by T. The daughter process 
corresponds to a Sibuya p.p. T(/i) with its every point given a multiplicity according to 
independent copies ofY^. 

Proof. In the case when the spectral measure a is supported by probability measures, 
representation (137]) becomes 

G*[h] =exp{ - f (1 - B{h),p) a a(dp)} Wh £ V(X), (38) 

JMi 

where Mx is the space of probability measures on X. In terms of the p.g.fl. (1271) of a 
Sibuya p.p., this reads 


[h] = exp { — / 1 — (1 — (1 — B(h), /i) Q )cr(d/i)} = 

./Mi 

= exp { ~ [ (l-Gr{ti)[B(h)])a(dn)} heV(X), (39) 

J Mi 

where T (pi) denotes a point process following the Sib(a,/i) distribution. Notice that, 
since by ([5]), B(-) is the p.g.f. of the distribution Y a Q , Gr^)[B(h)] is the p.g.fl. of a point 
process by (ED. □ 

As we have seen in (122]) . 2 r -stable processes are infinitely divisible. The latter can be 
divided into two classes: regular and singular depending on whether their KLM-measure 
is supported by the set of finite or infinite configurations (see, e.g., [5] Def.l0.2.VI]). 
Similarly to the proof of Theorem 29 in [7] on the decomposition of DaS processes, we 
can extend this result to Testable processes. 
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Theorem 5. An T-stable p.p. $ with a spectral measure a can be represented as the 
sum of two independent T-stable point processes: 

$ = 

where <f> r is regular and singular. <f> r is an T-stable p.p. with spectral measure being 
<j | Mi = cr(- nMj) and $5 is an T-stable p.p. with spectral measure ^|s\ Ml ■ 

The regular component <h r can be represented as a cluster p.p. with p.g.fl. given by 
G^ r [h] = exp { - f (l - G r ^)[B(h)])a\ Ml (dfj,)} Vh e V(X). 

J Mi 

On the contrary, the singular component is not a cluster p.p., and its p.g.fl. is given 
by (CUD (with a replaced with 0 "|g\ Ml there). 

4 General branching stability for point processes 

Stable distributions appear in various limiting schemes because decomposition of a sum 
into a proportion t and 1 —t of the summands inevitably leads to the limiting distribution 
satisfying ([I]). We have seen that this argument still works for point processes when 
the multiplication is replaced by a stochastic branching operation, the reason being 
associativity and distributivity with respect to the sum (superposition). One may ask: 
to which extent one can generalise this stochastic multiplication operation so that it still 
satisfies associativity and distributivity? The answer is given in this section: branching 
operations are, in this sense, the exhaustive generalisation. 

4.1 Markov branching processes on J\f 

Markov branching processes on Af (also called branching diffusions or branching particle 
systems) basically consist of a diffusion component and a branching component: each 
particle, independently of the others, moves according to a diffusion process and after 
an exponential time it branches. When a particle branches it is replaced by a random 
configuration of points (possibly empty, in which case the particle dies) depending on 
the location of the particle at the branching time (e.g. m. 0 or uni)- 

Alternatively branching particle systems can be defined as Markov processes satisfy¬ 
ing the branching property, as follows. 

Definition 7. A Markov branching process on Af is a time-homogeneous Markov process 
(f&t)t>o,ip<zN 011 (A/", B(Af)), where t denotes time and <p the starting configuration, such 
that its probability transition kernel Pt{<p, •) satisfies the branching property: 

Pt(<Pi + <P 2 , ■) = p t(<Pir) * Pt(<P 2 ,-), (40) 

for any t > 0 and ipi, p >2 in Af. 
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The branching property fl40j) can also be expressed in terms of p.g.fl.’s as follows: 


Gf[h] 


1 , if tp — 0, 

n xep G S t’[h], 


h e V(T), 


(41) 


where Gf and G s t x are the p.g.fl.’s of and respectively (see, e.g., [TJ Ch. 5.1] or p, 
Ch. 3]). Under some additional regularity assumption, every Markov branching process 
on A f (defined as above) can be expressed in terms of particles undergoing diffusion and 
branching see, for example, m 

In general not every starting configuration p e Af is allowed. In fact when f consists 
of an infinite number of particles, the diffusion component could move an infinite number 
of particles in a bounded set. Therefore in general one needs to consider only starting 
configuration f such that Gf[h] < oo for any t > 0 (for more details see, e.g., [ 8 | Ch. 5] 
or |TQj Ch. 1.8]). 


4.2 General branching operation for point processes 

Let • : (t, $) —» t • <f>, t G [0,1] be a stochastic operation acting on point processes on 
X or, more exactly, on their distributions. We assume • to act independently on each 
realisation of the point process, meaning that 

P(t«$ £ A) = [ P(fip e A)P*(dp), A e B(Af), t e (0,1], (Al) 
Jjt 

where P$ is the distribution of <f>. We require • to be associative and distributive with 
respect to superposition: for any t,ti,h £ ( 0 , 1 ] and <f>, <l>i, $2 independent p.p.’s on Af 

ti • (t 2 • $) = {tit 2 ) • $ = (A2) 

^•(^>1 + ^ 2 ) = (A3) 

where in (jA2jl and (jA3j) the different instances of the • operation are performed inde- 
pendently. Note that (A3) implies that t • 0 = 0 for all t e (0,1], where 0 is the empty 
configuration. 

Remark 3. Because of (lAljh • is uniquely defined by its actions on deterministic config¬ 
urations f G Af. In fact, given (1A1D . 4>, , 4 > 2 in HA 2 D and (1A3|1 can be replaced with 

deterministic point configurations ip, f\,fi £ Af. Note that, although p is deterministic, 
t • p is generally stochastic (as, for example, for the thinning operation). 

As we have seen in (120|1 and (12T| . the (local) branching operation o T operation satisfy 
m-m- The following results characterises stochastic operations satisfying (IA1D - (IA3I) 
in terms of Markov branching processes on Af. 

Definition 8. We call a stochastic operation • acting on p.p.’s on X a (general) branch¬ 
ing operation if there exist a Markov branching process (Tf )t>o,tpej\r on (A f, B(Af)), such 
that for any p.p. $ on A 

Tf = e“'«<I> t E [0, + 00 ). (42) 
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Proposition 1. A stochastic operation • satisfies ED-Ei if and only if it is a general 
branching operation. 

Proof. Necessity: Let • satisfy (IA1D - (IA3D . Let Pt(p>, •) denote the distribution of e~ l • ip. 
By putting $ = e~ tx • ip, ip E Af, in (ED we obtain 

P[e -i2 • (e~ tl •ip) £ A] = [ P t2 (ip,A)P tl (ip,d(p), A G B(Af), h,t 2 > 0. (43) 

J M 

Using the associativity of • (Assumption (lA2jl ) on the left-hand side of (1431) we obtain 
the Chapman-Kolmogorov equations 

p t 1+ t 2 (p,A) = [ P tl (ip,A)P t2 (ip,dip) A g B(Af), ti,t 2 > 0. (44) 

J M 

Therefore, by the Kolmogorov extension theorem there exists a Markov process x I’f on 
Af having transition kernel (Pt(<p, -)) t > 0 - Let p E Af\0 and Gf[-] be the p.g.fl. of Pt(g>, •) 
for t > 0. Using the definition of Pt(ip, •) and the distributivity of • (Assumption (|A3p ) 
we obtain 

G?M = G,-,. v [h] = g e , b «-, a [ft] = n G e -..Ufc] = nef-W, h e W 

xSLtp x€.tp 

From distributivity it also follows that e _i • 0 = 0 and therefore G®[h\ = 1 for any 
h G V(X). Therefore (l4Th is satisfied and Tf is a Markov branching process on Af. 

Sufficiency: Let (Tf ) t > 0 , <p G Af be a Markov branching process on Af with transition 
kernel P t (ip,-). Consider the operation • induced by (1421) . namely t • $ = v k $ logf . 
Assumption (IA1D follows by the construction: 

p(Tf g A) = [ P(vkf g A)P*(d<p), A g J3(Af), t > 0. 

J N 

Given <p G Af and ti,t 2 G (0,1], using flAlh and the Chapman-Kolmogorov equations 
(1441) we obtain 

P (U • (t 2 • <p) G A) = / P_ log tl ftp , A) P_ log t 2 (tp, dip) = 

JN 

— -P-iog(i 1 i 2 )( < ^) A) = P{{titf) • <p E A) A g B(Af), 
i.e. the associativity 0A2[) ) of • holds. 

Finally, let Gf[-] be the p.g.fl. of Tf for t > 0 and <p G Af. Given t G (0,1] and 
<Pi,<P 2 £ ff\ 0, using the independent branching property (HU), it follows that 

G t .(„ +M) [fe ]== n G H fe ]=n G ‘>’W n = 

x£<P2 

= Gy ogi [ft]Gy og ,[h]=G,. OT [h]G,., J [fc] heV(X). (45) 
The distributivity (IA3D of • follows from (1451) and Remark [3l □ 
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Example 3 (Diffusion). Let (X t ) t >o be a strong time-homogeneous Markov process on 
X , right continuous with left limits. Let be the Markov branching process 

on A f where every particle moves according to an independent copy of X t , without 
branching (see [U Sec. V.l] for a proof that this is indeed a Markov branching process 

on A/"). Denote by »d the associated branching operation, t $ = v h $ logt . Let Pt(x, •) 
be the distribution of X£, where x denotes the starting state, and Pth{x ) = E/i(Xf) = 
f x h(y)P t (x, dy ). Then, for any <p in A F we have 

a ~(z<p 

= J] Eh(X x _ logt ) = P-\o gt h{ X ) = G„[P_ logt /z] h G V(X). (46) 

x(E.<p 

Example 4 (Diffusion with thinning). Let X t be as in Example |4j Let (4/ t)t>o,<peAT 
be the Markov Branching process on A/", where every particle moves according to an 
independent copy of X t and after an exponentially Exp(l)-distributed time it dies (inde¬ 
pendently of the other particles). We denote by •dt the associated branching operation 
t m di t T = v L?io gt - Similarly to (l46jh it is easy to show that given a p.p. <f> we have 

Gfr dt $[h] = G$[l — t + t(P-\ ogt h)], h G V{X). 

This operation acts as the composition of the thinning operation o and the diffusion op¬ 
eration »d introduced in Example [3] regardless of the order in which these two operations 
are applied. For any p.p. <f> 

t »dt $ = t » d (t o $) = t o (t •d <f>). 

In fact, since 1 — t + t(P_ logt h) = P_ logt (l — t + th) for any h in V(X), it holds 

[h] GtoQ [P— log th\ G*ij>[l t T i(P_ logfh)] 

= G$[P_i og t(l — t + th)] = Gfr d $\ 1 — t + th\ = Gto(t» d ®)[h]- 

4.3 Stability for general branching operations 

Proposition H] shows that branching operations are the only operations on point processes 
satisfying assumptions flAll) - (1A3l) . Such assumptions, together with the continuity and 
subcriticality conditions below, lead to an appropriate definition of stability, as Propo¬ 
sition [2] below shows. 

Definition 9. We call a branching operation • continuous if 

t • p p for t f 1 for every p G A/", (A4) 

where =>■ stands for weak convergence (or equivalently for the convergence in Prokhorov 
metric). Moreover we say that • is subcritical if the associated Markov branching process 
on J\f, tyf, is subcritical, i.e. if ET^A) < 1 for every x G X and t > 0. 
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Proposition 2. Let $ be a p.p. on X with p.g.fl. G$[-] and • be a subcritical and 
continuous branching operation on X. Then the following conditions are equivalent: 

1. V n G N 3 c n G (0,1] such that given $( n ') independent copies of $ 

$^ Cn #(<I> (1) + ... + $ (ri) ); (47) 

2. V A > 0 3 t G (0,1] such that 

G*[h] = (GW[/i]) A ; 

3. 3 a > 0 such that VneN, given independent copies of & 

= (n _ “) • + ... + <h (n) ); (48) 

4- 3 a > 0 such that V t G [0,1] 

G*[h] = (G t .' S ,[h]y~°‘) (49) 

5. 3 a > 0 such that V t G [0,1] ; given (E^ 1 ) and d^ 2 ) independent copies of 

t 1/a • $ (1) + (1 - t) 1/a • $ (2) = $. (50) 

Proof. If $ = 0 then all the conditions are trivially satisfied. So we suppose $ ^ 0. 
4) => 2) =>■ 1) are obvious implications. So if one proves 1)=>4) then 1), 2) and 4) are 
equivalent. 

To show 1)=>4) n °t e that, given n G N, the coefficient c n satisfying (1471) is unique. 
In fact if c n and tc n both satisfy (1471) . with t G (0,1), by associativity it follows 

$ = ( tc n ) • (4> (1) + ... + 4> (n) ) = t • (c n • (4> (1) + ... + 4> (n) )) =t»$ 

and thus, because of subcriticality, t — 1. Using (1471) and the distributivity and asso¬ 
ciativity of • we obtain that, given m,n G N, 

4> = c n • (4 >( ' 1) + ... + $(”)) = 

= c n • (c m • (<h (1) + ... + $^) + ... + c m • ($("-i)"*+i + 

V ! \ 

— \CnCm) * 


... + $( nm ))) 
($(!) _|_ ... _|_ 


which implies that 

Cnm C n C m . 

Since we are considering the subcritical case, we have 

n > m =>■ c n < c m . 


(51) 

(52) 
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For every 1 < m < n < +cx), m,n G N define a function c : [1, +oo) D Q —* (0,1] by 
setting 



(53) 


The function c is well defined because of (15T]) and it takes values in (0,1] because of 
(]52]h Using associativity, distributivity and (l47jh 


n n 

( G sTT“ = (G e . (cm ., 4ll)+ +tW)) W)" = 

n 

= (G c „. ( tm + .„ +t W)(ft]) S = ((G 0 „.*[ft]) m ) ” = (G c „. s .[/ I ])" = (54) 

Therefore 

G$[h\ = (G c ^,^[h]) x Vx G [1, Too) nQ. (55) 

It follows from (J55J) and (155]) that c is a strictly decreasing function on [l,+oo) flQ. 
Therefore we can be extended to the whole [1, +oo) by putting 

c(x) := inf{c(y) : y G [1, x\ D Q}. 

From (1511) and (1521) . taking limits over rational numbers, if follows that c(xy) = 
c(x)c(y) for every x,y G [1, +oo). The only monotone functions c from [1, +oo) to (0,1] 
such that c(0) = 1 and c(xy) = c(x)c(y) for every x, y G [1, +oo) are c(x) = x r for some 
r G M. Since our function is decreasing then r < 0. Let a > 0 be such that r — — 1/a. 
Fix x G [1, +oo) and let {arjnGN C [1, +oo) D Q be such that x n l x as n —> Too, and 
therefore Xn 1 ^’ f x -1 ^ as n —y Too. Since • is left-continuous in the weak topology 
(assumption (1A4j) ) it holds that 

x~ 1/a • <F =>• x~ x l a • <F n —> Too, 

which implies 

G x -v«.M n -^ G x - 1/a .4h} V/ l GV(4 

From (1551) we have 

(G*[h]) 1/x = lim G c (a. n ).*[/i] = lim G -i/a [h] V h G V(A'), 

n—>+00 n—>+00 Xn 1 

and therefore we obtain (1551) as desired. 

4) => 3) => 1) are obvious implications and thus 3) is equivalent to 1), 2) and 4)- 
To show 4) => 5) take x, y G [1, Too). Then, because of 4), 


G,[h] = G {x+y) -^[hr y = G x _ UaM -„.Jh]* • G y _ H ^ yl/a Jh]y = 


= G 


( 4 s ) 




[h], (56) 


where is an independent copy of <f>. Then 5) follows since x, y G [1, Too) are arbitrary. 
5)=>3). (158]) can be obtained iterating (150]) n-1 times. □ 
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Section [4] shows that branching operations are the most general class of associative 
and distributive operations that can be used to study stability for point processes. There¬ 
fore the following definition generalises all the notions of discrete stability considered so 
far. 

Definition 10. Let <f> and • be as in Proposition [2j If (I5U]1 is satisfied we say that $ is 
strictly a-stable with respect to • or, simply, branching-stable. 

In this paper we do not provide a characterisation of stable point processes with 
respect to a general branching operation. Instead, next we consider some specific cases 
that point towards directions to obtain such a characterisation in full generality. The 
main idea is that, given a branching operation • acting on point processes, there is a 
corresponding branching operation 0 acting on random measures such that stable point 
processes with respect to • are Cox processes driven by stable random measures with 
respect to 0 . 

4.4 Stability with respect to thinning and diffusion 

Cox characterisation Recall that DoS point processes are Cox processes driven by 
StaS intensity measures, see Section 12.41 The main reason for this is that the thinned 
version of a Poisson p.p. with intensity measure p, 11^, is itself a Poisson p.p. with 

T> 

intensity measure tp, i.e. t o II M = fl^. The same holds for a Cox p.p. II^ driven by a 
random measure 

t°n 5 = n t€ . (57) 

For the thinning and diffusion operation •dt of Example [2 

G t . dt w IL [h] = Gnjl ~ t + tP-i og th] = exp{-(l - (1 - t + tP- logt h), p)} = 

= exp{-(LP_ logt (l - h),p)} = exp{-(l - h,tPf l t p)} = G UtP , [h], (58) 

-logt^ 

where P* iogt is the adjoint to the linear operator P-i Q gt- If we denote by Qdt the 
following operation 

Qdt. : (0, 1] x M —> M 

tQdt fi:=tPl logt fi, 

then (1581) implies t »dt II= U tQdtfl . Similarly for Cox processes, 

t *dt n e = n tQdt £, ( 60 ) 

where the operation Qdt acts on each realisation of £. The analogy between (157)1 and 
QgQD suggests the following result. 

Proposition 3. A point process is strictly a-stable with respect to »dt If and only if it 
is a Cox process 11^ with an intensity measure being strictly a-stable with respect to Qdt- 
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Proof. Sufficiency. Suppose £ is strictly cx-stable with respect to Qdt- Then using (1601) 
and the stability of £ 




-6V“, 


* (It 


tt// — tt/ I T\ n — TT — TT 


where the 1 If and 11" are independent copies of lb. Therefore Id? is strictly ct-stable 
with respect to •at- 

Necessity. Suppose $ is strictly ct-stable with respect to • at■ From (148]) we have that for 
any positive integer m 


G,[fc] = (a m - y .. dl 4h]) m = (g*[i 



h e V(X). (61) 


We need to show that G$[l — u], as a functional of u G BM + (X), is the Laplace 
functional of a random measure £ that is strictly ct-stable with respect to Qdt- Since 
1 — u may not take values in [0,1], the expression G$[l — u) may not be well defined. 
Thus we use (RID and define L[u) as 


L[u] = (g 4 [1 - nG 1/ °Pi.«„«l)’", u £ BM + (X), 


(62) 


noting that for some m big enough 1 — m 1// “PiogmM takes values in [0,1] and the right- 

a 

hand side of (162]) is well defined. Arguing as in the proof of Theorem 3.5 one can prove 
that L is the Laplace functional of a random measure £, Lg. Finally, £ is strictly n-stable 
with respect to Qdt. because for any u in BM + (X) 


L te d A u ) = E exp{ —(w, t Q> d t £)} = Eexp{-(u, LP* logt £>} = 

Eexp{ ftP— jog tU, £)} [fP_ i 0 g t u ), 

and, supposing u < f _1 , and thus (1 — tP_i ogt u) G V(ff), we have 

L^[tP_i ogt u] = G$[l - LP_i og Ai] = G^[l-t + tP_i ost (l-u)] = 

G t . dt<s> [l-u] = G^l-uf = L^uf. (63) 

Similar calculations also apply to the general definition of L[u] in (1621) . which includes 
the case u > t -1 . The fact that £ is strictly a-stablc with respect to ©<# follows from 
(1651) arguing, for example, as in (155|) . □ 


Levy characterisation and spectral decomposition Proposition [3] characterises 
stable p.p.’s with respect to •<& as Cox processes driven by stable random measures 
with respect to Qdt■ In this section we describe stable random measures with respect to 
Qdt in terms of homogeneous Levy measures (with respect to Qdt) and we show how to 
decompose such homogeneous Levy measures into a spectral and a radial component. 

Given A G B{M) and t G (0, 1] we define t Qdt A = {t Qdt /i : H G A}. The idea 
is to look for homogeneous Levy measures of order a with respect to Qdt, meaning that 
for any A in £>(A4) 

A(f Q dt A) = t~ a A(A) VtG(0,1], (64) 


23 


Let us consider Laplace functionals of the form 


L[h] — exp \ — f (1 - e - ^) A(dp)\, h G BM+(X), (65) 

l Jm\{o} J 

where A is a Radon measure on A4\{0} such that 

f (l — e - ^)) A (dp) < oo, (66) 

J M\{ 0} 

for any h in BM + (X), and (IMj) holds for any A in B(A4). Arguing as in the proof of 
Theorem 2 of p], it can be seen that (j65]) defines the Laplace functional of a random 
measure, say £. Then, defining L^[h\ as in (l65lb from (1641) it follows that 

Lt Qd Ah] = exp (- [ (l - e~^) t~ a A(dp)\ = h e BM + (X), 

l dx\{o} J 

which means that £ is ct-stable with respect to Q d t by an argument analogous to the 
one in (1561) . 

We now show how to decompose Levy measures satisfying (1641) in a radial component 
(uniquely determined by a) and a spectral component. Such a spectral decomposition 
depends on the operation Q dt and thus it is not the one used in Section [ZU for thinning- 
stable point processes. For simplicity we restrict ourselves to the case where Ai is the 
space of finite measures on X = M” for some n and the the diffusion process Pt is a 
Brownian motion, meaning that given fi in A4 and t in (0,1], the measure t Q dt /i is 

tOdtB = tu t *fjL, 

where * denotes the convolution of measures and u t , for t in (0,1), has the following 
density with respect to the Lebesgue measure 

du t 1 f \x\ 2 \ 

d? = /,w = (2n^jf exp \=2toii|’ 

while for t — 0, v t equals So, with 0 being the origin of M n . 

We now show that M. \{0} can be decomposed as S x (0,1], for the following S 

S := {/i G M. \{0} : $(t, p) G (0,1) x M. \{0} such that t Qdt P = p}- 


Note that § G £>( M) because S = U4 g (o,i)nQ^0dt (A / 1\{0}) and tQdt(-M \{0}) G B(M) 
for any t in (0,1) DQ. In fact if A G B^AA), then also tQ d tA G B[M) because the map 
p —> t Q dt p is injective (see Proposition |J] below) and therefore the image of a Borel set 
is still Borel (see, e.g., Section 15.A of (IS]). 
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Proposition 4. The map 


(0,1] x -M\{0} 
(t,p o) | —> t Q dt p 0 , 


(67) 


is a bijection. 

Proof. Injectivity: First note that given t G (0,1) and p, p G M. such that u t *p = v t *Pi 
then p — p. This follows, for example, by noting that, for any t in (0,1], the moment 
generating function of u t is strictly positive on M n . Then, given ti,t 2 G (0,1] and 
Poi Po G S 

h Qdt Po = (tit 2 ) Qdt Po t\(v tl * p 0 ) = (tit 2 )(v tlt2 * po) 

"t! * Po = v tl * (v t2 * (t 2 p 0 )) Po = v t2 * (t 2 p 0 ) <*t 2 = 1 and p 0 = Po, 

which means that the map (t, p 0 ) i-> t Q dt p 0 is injective. 

Surjectivity: Let p G A4 \ {0}. Without loss of generality we can consider p to be a 
probability measure, otherwise consider Define 

1^ := {t G (0,1] : 3p G M. \ {0}, t G (0,1] such that t Qdt. p = p}- 

Note that 1 G because 1 Q dt p = p and, thanks to the associativity of Q dt , is an 
interval. We define to inf 7 /t and we prove to > 0. Suppose to = 0. Then for any 
£ G (0,1] there is p*D g M \ {0} such that p = eu £ * p( £ \ Setting p( £ ) = ep (D we have 
p = u £ * p^ £ \ Since both p and p& by associativity of Qdt are obtained by convolution 
with some v t , they both admit bounded density functions, say, g and g( £ \ respectively. 
Moreover ||p^||i = 1 because 

1 = Nil = ll/s *P (£) lli = I [ fs(x-y)g (£ \y)dydx = [ g {e \y)dy = ||p (e) ||i, 

J R n JR n JR n 

were we used the fact that both f £ and g ^ are positive, Fubini’s theorem (which holds 
for positive functions) and ||i/ e ||i = 1 and. Then 

IIpIIoo = Wfe* g (£) \\oo < ||/e||oo||p (£) ||l = ||/e||oo ^ 0 US £ ->■ 0 . 

Therefore, by contradiction, t 0 > 0. Given t G (to, 1] let p® G M. \ {0} such that 
u t * pQ = p. We now prove that {p®} as t ft 0 is Cauchy in the Prokhorov metric. Let 
£ > 0 be fixed and 5 = 5(e) > 0 to be fixed later. Consider fi,t 2 such that f 0 < t\ < 
t 2 < to + 5 < 1. We have p ^ = v tl /t 2 * p^ 1 / As 5 —> 0 we have v tl / t2 (W. n \ B e ( 0)) —> 0. 
Thus we can choose 5 such that / <2 (M n \ B e ( 0)) < e. Therefore we have 

p {t2 \A) = u tl/ta * p [tl \A) = f u tl/t2 (A - y) p (tl \dy ) = 

J R" 

= [ v tl/t2 (A-y) p {tl \dy) + [ v tl/ta (A - y) p {tl \dy) < e + p {tl) (A £ ), 
Jm n \A e J A e 
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where A — y is defined as {x G X : x + y £ A} and A £ = {x G X : B e {x) D A ^ 0}. 
Thus there exists a probability measure G A1 \ {0} such that fi^ =4- fi^ as 111 0 
implying u t * =»■ zy 0 * /T to \ Therefore i/ to * = //. Finally we note that G § 

because of the definition of to- □ 

Thanks to Proposition [H for any // in A4 \ {0} there is one and only one couple 
(t, /io) G (0,1] x § such that t Qdt l 1 o = AC meaning that A4 \ {0} can be decomposed as 

M \ {0} = (0,1] x S. (68) 

Any measure A satisfying (jMj) can then be represented as A = 9 a ®a, where 9 a ((a, 6]) = 
( b~ a — a~ a ) for any (a, b] C (0,1] and a (A) = A ((0,1] x A) for any A G £?(§). In fact 
for any (a, b] C (0,1] and A G B(S) 

A ((a, b] x A) = A((0, b] x A) - A((0, a] x A) # 

= 6 _a A((0,1] x A) - a _Q A((0,1] x A) = 6 a ({a, 6]) <r{A). 

Since 9 a is fixed by a, there is a one-to-one correspondence between Levy measures 
satisfying (164|) and spectral measures a on §. Thus a Cox point process driven by the 
parameter measure £ is strictly a-stable with respect to the diffusion-thinning operation, 
if £ is ©^-stable, which, in turn, can be obtained by choosing an arbitrary spectral 
measure a on § satisfying 

fi(B) a a(dii) < oo 

for all any compact subsets 5cK“. And then taking £ with the Laplace functional 

L^[u] = exp j— J~( u j n) a &(dii )|, u G BM + (A). 

As an example, consider a finite measure cr on M n and a its push-forward under 
the map x G X i —y S x G A4. As shown above, the homogeneous Levy measure on 
A4 = (0,1] x S having a as spectral measure is A = 9 a © a. Therefore A is supported 
by the following subset of AT 

y := {v t { ■ - m) : t G (0,1], m G M n } C M. 

4.5 Stability on Z + and M + 

In all the examples considered so far: thinning stability, ^-stability and an example 
of a general branching stability, the operation in Section 14.41 the corresponding 
stable point processes were Cox processes driven by a parameter measure which is itself 
stable with respect to a corresponding operation on measures. Unlike its counterpart 
operation on point processes, this operation was not stochastic, meaning that its result 
on a deterministic measure is also deterministic. For the case of thinning stability, this 
was an operation of ordinary multiplication, for a point-processes branching operation 



26 



•dt, this was the operation Qdt- Nevertheless this does not need to be the case in general: 
in this section we consider point-processes branching operations • whose corresponding 
measure branching operation © is also stochastic. 

We consider the case of a trivial phase space X consisting of one point, or in other 
words the case of random variables taking values in Z + and M + . For Z + , since the phase 
space consists of one point, the general branching stability corresponds to the Testability 
described in Section 12.31 We introduce the notion of branching stability in M + using 
the theory of continuous-state branching processes (CB-processes, see [T7j|). We show 
that, at least in the cases we consider, branching stable (or Testable) discrete random 
variables are Cox processes driven by a branching stable continuous random variable. 
Finally we show how to use quasi-stationary distributions to construct branching stable 
continuous random variables. 

We argue that the theory of superprocesses (e.g. m) should be relevant to extend 
the ideas presented in this section to general branching stable point processes. 

Continuous-state branching processes Continuous-state branching processes were 
first considered in [13j and na, and can be thought of as an analogue of continuous time 
branching processes on Z + on a continuous space M + . 

Definition 11. A continuous-state branching process (CB-process) is a Markov pro¬ 
cess (Zf ) X)t >o on M + , where t denotes time and x the starting state, with transition 
probabilities (P t x ) tx>0 satisfying the following branching property : 

pp" = n * p; j , (es) 

for any t , x, y > 0, where * denotes convolution. 

A useful tool to study CB-processes is the spatial Laplace transform V t , defined by 

xV t (z) = -log / e~ zy P t x (dy ) z> 0, (70) 

J R+ 

for t > 0 and x > 0. The value of x in (1701) is irrelevant because of the branching 
property, it could be simply set to 1. Using the Chapman-Kolmogorov equations it 
follows from (1701) that (Ut)teR+ is a composition semigroup 

v t (v s {z)) = v t+s (z ), s,t,z > o. (cr) 

Similarly to the discrete case in Section 12711 we focus in the subcritical case, E^ 1 ] < 1, 
and we assume regularity conditions analogous to (1C2D - (1C4D . More specifically, rescaling 
the time by a constant factor if necessary, we may assume that 

E[Zl] = e~ l , (C20 

lini Vt(z) = Vq(z) = z, (C3 r ) 

lim V t (z) = 0. (C40 

t—> oo 
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Example 5. A well known example of CB-process is the diffusion process with Kol¬ 
mogorov backward equations given by 

du du bx d 2 u 

dt aX dx 2 dx 2 ' 

The spatial Laplace transform of the corresponding CB-process is 


Vt(z) 


z exp (at) 


1—(1—exp(at))!f ’ 


1 +t 


bz ? 


if a ^ 0, 
if a = 0. 


(71) 


The sub-critical case corresponds to a < 0. Rescaling time to satisfy (1C2'|) corresponds 
to setting a = — 1. 

The limiting conditional distribution (or Yaglom distribution) of ) x ,i>o is the weak 
limit of (Zf |Zf > 0), when t —> +oo. Such limit does not depend on x (e.g. [16j Th. 3.1] 
or [18, Th. 4.3]) and we denote by the corresponding random variable and by Lz x 
its Laplace transform. The Yaglom distribution is also a quasi-stationary distribution, 
meaning that > 0) = Z ^ (e.g. [HJJ Th. 3.1]). Given (1C2 / j) and the quasi- 

stationarity of Z^ it follows that 

L Zao {V t (z)) = l-e- t + e- t L Zoa (z), s,z> 0. (72) 


V-stability and Cox characterisation of Testable random variables Let (Z?) x,t> 0 
be a CB-process with spatial Laplace transform V = satisfying assumptions 

dcrj-tfcl} of the previous section. Dehne a corresponding stochastic operation acting 
on random variables on M+ as follows: 


t&vt = Z[ logt 0<z<l, (73) 

where ^ is an M + -valued random variable and Zf is the CB-process with random start¬ 
ing state £. Similarly to Proposition [Q from the Markov and branching properties of 
(Z*) Xtt > 0 it follows that ©y is associative and distributive with respect to the usual sum. 

The notion of V-stability for continuous random variables is analogous to the notion 
of ^-'-stability for discrete frameworks: 

Definition 12. A M + -valued random variable X (or its distribution) is V-stable with 
exponent a if 

t 1/a ©y X' + (1 - t) 1/a ©y X" = X 0 < t < 1, (74) 

where X' and X" are independent copies of A". 

In terms of Laplace transform L, the definition of ©y in (1751) can be written as 
L t Q v x(z) = Lx(V-\ og t(z)). Thus, arguing as in Proposition [2j (174|) is equivalent to 

L x (V- logt (z)) = L x (z) ta 0<t<l. (75) 
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Suppose we have a continuous-time branching process on Z + with p.g.f.’s T = ( F t ) t > 0 
and a CB-process with spatial Laplace transform V = (V t )t>o such that 

F t {z) = 1 - V t (l - z) 0 < 2 < 1. (76) 

The relation between the discrete and continuous stochastic operations, o T and ©y, is 
that Cox random variables driven by a V-stablc random intensity are J-'-stable. 

Proposition 5. Let T = (F t ) t >o and V = ( V t )t>o satisfy (1C1 jl ~ (1C4[) and (ICl'jbflCVj) 
respectively, and let (176)1 be satisfied. Let f be a V-stable random variable with exponent 
a and let X be a Cox random variable driven by f, meaning that X|£ ~ Po(f). Then 
X is F-stable with exponent a. 

Proof. The pg.f. of X is given by Gx(z) = L^( 1 — z), see e.g. (fT2j) . Therefore 

G x (F. Iosl (z)) = L ( ( 1 - (F_ logi ( 2 ))) 1* L ( (V- , ogl (l - z)) = L g ( 1 - zf = G x (z)‘°, 

which implies that A" is Testable with exponent a. □ 

Examples of discrete and continuous operations, o T and ©y, with T and V coupled by 
(I77)j) are thinning and multiplication, as well as the birth and death process of Example 0 
and the CB-process of Example [5] (with b = 1). Also the operations •dt and Qdt of 
Section l4~4l in a point process and random measures framework, satisfy fl76|h A natural 
question is whether for any continuous time branching process on Z + there is a CB- 
process such that (177)1) is satisfied and viceversa. Note that, given (176|) . (F t ) t > 0 is a 
composition semigroup if and only if (V t ) t >o is. Indeed, 

FfiFfiz)) = l-V t {l-(l-V,(l-z))) = l-V t {V 8 (l-z)) = l-V t+s (l-z) = F t+s (z), 

and similarly for V t . Therefore one would only need to prove that if V t is the spatial 
Laplace transform of a random variable on M + then F t defined by (176|) is the p.g.f. of a 
random variable on Z + or viceversa. 

Finally, note that Proposition[5]suggests that, at least in some cases, A'-stable random 
variables on Z + are Cox processes driven by V-stable random variables on M + . It is 
therefore natural to ask whether we can characterise V-stable random variables. 
Equations (1721) and (175|) imply that 

tol 

tQvZ^^Y.Zoo 16(0,1], (77) 

i =0 

where t o 1 is, by the definition of thinning, a binomial Bin(l,£) random variable (in¬ 
dependent of Zoo). Therefore the Yaglom distribution allows us to pass from Qy to 
thinning and use such a property to construct V-stable random variables from DaS 
random variables. 

Proposition 6. Let X be a DoS random variable on Z + (see (USD ), and € = YfiLi z ™> 
where zQ , Z& , ■ ■ ■ are i.i.d. copies of the Yaglom distribution Z^. Then f is V-stable 
with exponent a. 
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Proof. Given its definition, the Laplace transform of £ is given by L^(z) = G x (L Zoo {z)). 
Therefore 


mi 


L^(V. logt (z)) - G x (L Zoo (V_i ogt (z))) - G x (1 ~t + tL Zoo (z)) - 

cm 


G toX (L Zaa (z)) ^ G x (L Zoo (z)) = L&Y , 


which implies that £ if V-stable. □ 

5 Discussion 

In this paper we have studied discrete stability with respect to the most general branch¬ 
ing operation on counting measures which unifies all notions considered so far: discrete 
stable and Testable integer random variables, thinning-stable and J-’-stable point pro¬ 
cesses characterised above. We considered in detail an important example of thinning- 
diffusion branching stable point processes and established the corresponding spectral 
representation of their laws. We demonstrate that branching stability of integer random 
variables may be associated with a stability with respect to a stochastic operation of 
continuous branching on the positive real line and we conjecture that this association 
may still be true in general for point processes and its continuous counterpart, random 
measures. A full characterisation of the branching-stable point processes, as well as of 
the associated stable random measures is yet to be established. 
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